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The influence of substitutional disorder on the magnetic properties of disordered Heisenberg binary 
spin systems with long-range exchange integrals is studied. The equation of motion for the magnon 
Green's function which is decoupled by the Tyablikov approximation is solved in the Blackman- 
Esterling-Berk(BEB) coherent potential approximation(CPA) framework, where the environmental 
disorder term is treated by virtual crystal approximation. The long-range exchange integrals in- 
clude a power-law decaying and an oscillating Ruderman-Kittel-Kasuya-Yosida(RKKY) exchange 
interaction. The resulting spectral density, which is calculated by CPA self-consistent equation, is 
then used to estimate the magnetization and Curie temperature. The results show, in the case of 
the three-dimensional simple cubic systems, a strong influence of ferromagnetic long-range exchange 
integrals on the magnetization and Curie temperature of the systems, which is obviously different 
from the calculation of short-range interaction. 
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I. INTRODUCTION 

The effect of disorder has been shown to be crucial for 
many essential properties of magnetic materials. How- 
ever, investigations involving statistical disorder are al- 
ways faced with lots of technical difficulties, especially in 
a parameter regime that is not accessible to perturbation 
theorji The introduction and application of the coher- 
ent potential approximation(CPA), which is initially de- 
veloped by SovenS and Taylor—, has increased greatly our 
understanding of disordered systems, particularly of ran- 
dom, substitutional disordered alloys. In the CPA, one 
assumes that the real, disordered material is replaced by 
a translationally invariant effective medium. This is fixed 
by imposing the self-consistent condition that the scat- 
tering off of a real atom embedded in the medium van- 
ishes on the average. The CPA provides reliable results 
in various situations studying excitations like electrons, 
phonons, magnons, eto^&. In recent years, CPA com- 
bined with some ab initio method has been shown to 
be crucial tool for the understanding of diluted magnetic 
semiconductors(DMS)i&. 

For spin systems with substitutional disorder, the 
problem more generally considered is a binary system^ 
in which the magnetic atoms of each component (A or 
B) with spin Sa or Sb are randomly distributed over the 

sites of a regular lattice with concentrations ca and 
cb (ca + cb = 1), respectively. In the present work, the 
interaction between spins is described by the Heisenberg- 
type Hamiltonian in which exchange constants are : 
A, A' = A or B according to the two interacting spins at 
site i and j. Using Tyablikov decoupling approximation 
to the retarded magnon Green's function of a Heisenberg- 
type Hamiltonian, we can achieve the equation of motion 
for Green's function involving three kinds of disorder that 
depend explicitly on site spin (diagonal disorder), the ex- 
change parameters (off-diagonal disorder) and an envi- 



ronmental term given by the static field induced at one 
site by the presence of neighboring spins (environmen- 
tal disorder). We will deal with diagonal disorder and 
off-diagonal disorder terms in the BEB-CPA framework. 
The BEB theory by Blackman, Esterling and BerkiS ex- 
tends the CPA which studies systems only with diago- 
nal disorder to random off diagonal matrix elements of 
the Hamiltonian (off-diagonal disorder) by doubling the 
Hilbert space (for a binary alloy) and adding the atomic 
sort to the matrix indicesii^. The environmental disor- 
der term is treated in the virtual crystal approximation 
in the same manner as done by Theumann*(to be re- 
ferred to as I), where an atom A(B) will be acted upon 
by a mean field so that full randomness will be preserved 
in the locator expression of CPA. Furthermore, consider- 
ing the Tyablikov decoupling approximation for the dis- 
ordered Green's function, magnetization of the different 
constituents need to be determined self-consistently for 
a given temperature because the locator and the local 
magnon spectral density are temperature dependent. In 
general, the analytical form of the local magnon spectral 
density can be got only for some simple cases such as the 
exchange integrals restricted to only nearest neighbors^. 
Analytically, the calculations where Jy are long ranged 
is not an easy task. However, long-range interaction are 
always of interest in different fields of physics because 
they can give rise to a variety of unusual macroscopic be- 
havior. For example, metallic magnetic materials are of 
itinerant type, which means that the exchange integrals 
between different localized magnetic ions are long range 
and driven by the polarization of the conduction carrier 
gas, i.e. the Rudcrman-Kittcl-Kasuya-Yosida (RKKY) 
type exchange interaction^*. 

In the present work, we calculate magnetization and 
Curie temperature of different components with fer- 
romagnetic short-range, power-law decaying long-range 
and oscillating RKKY long-range exchange interaction. 
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It should be mentioned here that the idea of using the vir- 
tual crystal approximation and BEB-CPA to deal with 
three kinds of disorder (diagonal, off-diagonal and envi- 
ronmental) of a binary spin system is same as paper I. 
However, our approach to solve the self-consistent CPA 
equations, which differs from paper I, follows the BEB- 
CPA method. Moreover, we present the method to cal- 
culate the temperature dependence of magnetization and 
the Curie temperature of systems by combining the meth- 
ods of paper I and the RPA-CPA theory^ done by G. 
Bouzerar et. al. Significantly, we can calculate any 
range exchange integrals (including short-range and long- 
range) . 

The article is organized as follows: In section we 
deal with different types of disorder and derive the config- 
urational averaged Green's functions, magnetization and 
Curie temperature for systems with ferromagnetic short- 
range, ferromagnetic long-range and RKKY-type long- 
range exchange interactions. The numerical calculations 
of the self-consistent equations are discussed in lllll for ex- 
change integrals of a long-range power-law decay | Jr\ ~ 
R~ a and a RKKY-type \J R \ ~ (RcosR - sini?)/i? 4 . In 
section llVl we will conclude the article with a summary. 



II. THE MODEL 



To study the magnetic properties of a material A CA B Cfl 
carrying a localized magnetic moment at each lattice site, 
we assume the dynamics of these spins to be reasonably 
described by the isotropic Heisenberg Hamiltonian 



H 



S 7 * • So 



i 



(1) 



where subscripts i and j refer to the occupied lattice 
sites, and Si = (Sf , Sf , Sf) is the spin operator of the 
localized magnetic moment at lattice site i with lattice 
vector R,j. The Hamiltonian contains a Zeeman coupling 
of the spins, where the external magnetic field is B — 
(0, 0,B). The effective Heisenberg exchange parameters 
Jij obey Jij — J(|Ri — Rj|) and Ju — restricting the 
theory to uniform media and disregarding more general 
case. Firstly, we generalize the Hamiltonian such that 
Jij is random: 



J, 



jAA 



if sites i, j are of type A 
if they are of type B 
= jf: B = jf, A if one is A site and the other B. 
Now introduce occupation indices Xi and yi, such that 
Xi = l,yi = if i is an A site, 
Xi = 0, yi = 1 if i is an B site. 
Some examples of their properties are: 

Xiy% — 0, — Xi, 

(Xi) c = c A , (yi) c = c B , 

where the bracket (...) means a configurational average. 
Then, we get 



J, 



, j — XiJ^ Xj + Xi Jjj yj + yi Xj + yi yj (2) 



Introducing the retarded magnon Green's function 

G l3 {E) = ((St-Sj))T (3) 

built up by the step operators Sf = Sf±iSf, its equation 
of motion reads 

(E - gjn B B) G l3 (E) = 2h 2 5 tJ (Sf) - 2h^^ m 

m 

x {{{SfS z m ; SJ)) T - ((S+Sf; S-)) «*) . (4) 

Making the Tyablikov approximation, which consists in 
decoupling the higher Green's function on the rhs. of J3J, 
one obtains after rearrangement: 

[E - gjH B B - 2h^2 Jim (S^)]G l3 (E) 

m 

= 2h 2 (Sf) % - 2h (Sf) JimG mj . (5) 

m 

Substituting Ji m by equation (J2J, we use the various 
combinations {x,y} to multiple the equation j£| in BEB 
manner. Then, the Green's functions' equation of motion 
can be expressed as, 



Si (My 



(6) 



if we introduce the 2x2 matrices 
/ qAA qAB 



r^BA s~f 



^ B * ~ [yiGiiXj yi G m ) ' ^ 



jAA rAB 

tBA jBB 

ij ij 



and 



Si = 



(xi9 A 


o x 


)-( 


\ o 


Vi9 B y 





Vi-Z* 



where a\ = (Sf ) /<jq (A = A or B, er = ca(S a ) + 
c B (S B )), w = (E- gjfi B B)/2h<7 and 

"A = X)( J ™ W X ™ + 4n (Sfn) Vm)/^ , (8a) 
m 

"b = ]T( J ™ ( S m) x m + J BB (Sf n ) y m )/a . (8b) 

m 

Now we find, as in BEB, that the matrix propagator 
&ij satisfies an equation for a problem with only local 
disorder, the matrix being independent of the random 
numbers. However, the so called 'local potential' u>a and 
lob are functions of random variables. Using the virtual 
crystal approximation as done by paper I, the equation 
© can be written as 

WA = C A OTA(Y J J «n) + CB<T B (Y, J fn?) , M 

m m 

lob = c A a A C£j^) + cWE J - ) • ( 9b ) 
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Considering, for any site i, J^m^irn — constant in an 
exchange-integral and lattice-structure given spin sys- 
tem, the 'local potential' uja and u>b, which describe local 
site occupation, are constants. In other words, the ex- 
change interaction of neighbor spins of site i leads to a 
static field which acts on site i. The static field is of 
independence on the location of site i as the external 
magnetic field B does. Thus, in our approach, the ex- 
change interaction of neighbor spins of site i is taken as 
an additional static field. The difference between uja and 
lob root in whether an atom of type A or of type B occu- 
pies site i. In our model, the random A- type or B-type 
occupation of any site decide the 'local potential' uja or 
lub, respectively. From here, our calculation will follow 
the BEB-CPA methods and differ from paper I. 

If we refer to three specific cj-dependent parameters as 
a matrix form in BEB-CPA manner 



Un = 



Ui u 3 
u 3 u 2 



the CPA self-consistent condition is expressed asAi 

7=(G 4l ) c = 7V- 1 ^G k , (io) 

k 

where Gk = [7 _1 + Uo + Jk] _1 and 

_ f.J AA (k) J AB (k) 
Jk ~ \j BA (k) J BB (k) 

where J AA (k) is the Fourier transform of the parame- 
ters of exchange interaction J AA (R). The renormalized 
locator can be written as 



7 



7 







where 



cao-a 



U) — LUA ~ O ' aV\ ' 

B _ c b <tb 

U) —U>B — ObU-2 



(11a) 
(lib) 



In general, the introduced 2x2 matrix Uo describes ef- 
fective medium propertiesSaii. From the denominator of 
the renormalized locators (|11|) . we find Uo has the same 
dimensions of exchange parameters J. So, according to 
BEB-CPA theory, the effective medium parameters Uo 
present the appropriate sum of the exchange integrals 
and are called the renormalized interactor. 

Once we get two renormalized locators j A and j B from 
CPA self-consistent equations, the magnon spectral func- 
tion can be expressed as 



1 Im7 A (w) 

7T 



(12) 



Furthermore, magnetization reads 

. (S x - $ A )(1 + 4>a) 2 ^ +1 + (1 + S x + $ A )$f -+ 1 



(Si) 



(1 + $a) 2Sa+1 



2S> + 1 
A 



(13) 



where A = A or B and the average magnon number $a 
can be calculated by 



$A = 



du> 



e 2ha ^/k B T _ 1 



(14) 



Now, for a given temperature and good starting value 
of {S%), we can get the renormalized locators 7 and 
the magnon spectral functoin D\(uj) from CPA self- 
consistent equations (|10|l . Furthermore, by using the 
equation Q13[l. we can calculate new values of magne- 
tization (Sf), which are re-inserted in a\ and 7 A . 

Considering (S%) — > in the limit T — ► Tq, we can get 
the Curie temperature 



2 2Sr S X (S X + 1) 

k B T c = -n L CA Fx 



where 



(15) 



(16) 



To study the influence of the range of a ferromagnetic 
exchange interaction on magnetization, we need to cal- 
culate J AA (k) and Ylm^im with ferromagnetic decaying 
long-range 



J XX '(R) 




(i?/a)~ 4 for R<R cut -off, 
otherwise. 



(17) 



Here, -R C ut-off is the maximum distance out of which the 
spin-spin interaction is not taken into account in the cal- 
culations, a is the lattice constant and Jo the nearest- 
neighbor interaction strength. For ii cu t-off = a, we get 
the case of nearest- neighbor interaction J = Jo, other- 
wise J = 0. At the same time, in the numerical studies, 
we also use RKKY-type long-range exchange interaction 



J AA '(i?) = Jo X C~) 



sin(-) - (-)cos(-) 
a a a 



(18) 



Considering J (i?) depending on the shell R, we have 



£j ^=]>>J^(^) 



(19) 



where Y]j corresponds to a summation over the I-th shell 
with a distance Rj from a given site and zj is the total 
number of spin in the shell. We also get a shell expression 



J AA '(k) =J2 J^'iRl^M^m), (20) 



where the sum runs over each site located in the I-th 
shell. 
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FIG. 1: Magnetization (S z ) as function of temperature T for 
nearest-neighbor interaction (Rcut-off ~ a for power-law de- 
caying long-range exchange integrals) on a simple cubic lattice 
for various concentrations ca- The parameters are J AA = 0.2, 



0.1, J<f 



0.15, S A = 2 and Sb = 1- 



III. NUMERICAL STUDIES 

In this section, we will mainly consider the influence of 
long-range exchange integrals on Curie temperature of a 
three-dimension disordered binary simple cubic systems 
with different exchange constants. For simplicity, we only 
consider the case of a zero external field and the lattice 
constant equals to 1. 

The magnetization of two different spin components 
as a function of temperature is shown in Fig. ^for sev ~ 
eral different concentrations ca- The figure shows that 
each spin component's magnetization equals to their spin 
quantum number at T = and will decrease monotoni- 
cally with increasing temperature. For any concentration 
ca, (Sa) an d (S B ) will decrease to zero at the same tem- 
perature. It is because the exchange integrals between 
two components is ferromagnetic and the ferromagnetic 
ordering will disappear at the Curie temperature. 

In Fig. [3 we show the Curie temperature as a func- 
tion of the maximum length R cut - ff of exchange inter- 
action, i.e. the spin-spin interaction will not be taken 
into account out of R cut - ff. For power-law decay- 
ing long-range exchange integrals, it can be found in 
Fig. [5] that the Curie temperature increases monotoni- 
cally with increasing R cu t-off and has a saturation value 
at Rcut-off ~ > oo for any concentration. For RKKY- 
type long-range exchange integrals, the Curie tempera- 
ture firstly increases with increasing R cut - ff and then 
presents an oscillating behavior around a fixed value, 
which corresponds to R cu t- ff — > oo. Both calculated 
results in Fig. [21 show that long-range exchange integrals 
will lead to changes of the Curie temperature and mag- 
netic properties of spin systems compared to the short- 
range case. The results also show that the long range of 
the interaction has to be taken into account explicitly in a 
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FIG. 2: Curie temperature Tc as function of the effective 
radius Rcut-off /a (a: lattice constant) for (a) power-law de- 
caying and (b) RKKY-type long-range exchange integrals on 
a sc lattice. The parameters are Jq A = 0.2, J AB = 0.1, 
Jo B = 0.15, Sa = 2 and Sb = 3. 



respective model calculation. To cut at a too early stage, 
for mathematic simplicity, may lead to rather misleading 
results. 



Furthermore, considering various ratios of ca to cb 
that affect Tc, the Curie temperature as a function of 
ca for nearest-neighbor and infinite long-range exchange 
integrals (R cut -off —> oo) is shown in Fig. [31 At the 
same time, the figures show the influence of three differ- 
ent strengths Jq B on the Curie temperature. Note that 
pure A (or B) corresponds to ca = 1 (or ca = 0). Com- 
pared with other methods^, the Curie temperature in 
the left column of Fig. [21 which is the result of nearest- 
neighbor exchange parameters, present similar tendency 
for different values of J$ B . The spin Sa, Sb and ex- 
change integrals J AA (A, X' — A or B) of the binary spin 
system will influence magnetization and Curie tempera- 
ture corporately. Comparing the left column to the right 
column of the Fig. [3 we find that the Curie tempera- 
ture of long-range exchange integrals is hugely different 
from the Curie temperature of nearest-neighbor exchange 
integrals. 
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FIG. 3: Curie temperature Tc as function of the concentra- 
tion ca for power-law decaying and RKKY-type long-range 
exchange integrals on a sc lattice, where nn means nearest- 
neighbor exchange parameters and inf means R C ut-off — » oo. 
Note the different scales for Tc at the left and the right col- 
umn. The parameters are Jo A = 0.2, J SS = 0.15, Sa = 2 
and Sb = 3. 



ing the Tyablikov approximation to equation of motion 
for the magnon Green's function, we investigate substi- 
tutional random spin systems in BEB-CPA framework. 
The resulting theory is then solved numerically in a self- 
consistent way for a simple cubic systems. The obtained 
spectral densities are then used to calculate spontaneous 
magnetization and estimate the Curie temperature of the 
corresponding system. 

For ferromagnetic power-law decaying and RKKY- 
type long-range exchange interactions, the calculations 
show long-range exchange integrals will lead to changes 
of macroscopic properties such as the Curie temperature 
of spin systems. The Curie temperature will oscillate 
for RKKY-type long-range exchange interaction and in- 
crease monotonically for power-law decaying interaction 
with the effective range R cu t-off of exchange integrals. It 
will tend to a saturation value for R cut - ff —* oo. These 
results imply that the effective exchange interaction, nu- 
merically taken into account, must be very long-ranged. 
The influence of long-range interaction on magnetization 
and Curie temperature should carefully be considered in 
the theoretical model. 
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